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We study the detection of weak coherent forces by means of an optomechanical device formed 
by a highly reflecting isolated mirror shined by an intense and highly monochromatic laser field. 
Radiation pressure excites a vibrational mode of the mirror, inducing sidebands of the incident field, 
which are then measured by heterodyne detection. We determine the sensitivity of such a scheme 
and show that the use of an entangled input state of the two sideband modes improves the detection, 
even in the presence of damping and noise acting on the mechanical mode. 
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I. INTRODUCTION 

Optomechanical systems play a crucial role in a variety 
of precision measurement like gravitational wave detec- 
tion pj and atomic force microscopes 0]. These systems 
are based on the interaction between a movable mirror, 
the probe experiencing the force to be measured, and a ra- 
diation field, the meter reading out the mirror's position, 
which is due to the radiation pressure force acting on the 
mirror. The mechanical force exerts a momentum and 
position shift of a given vibrational mode of the mirror, 
which in turn induces a phase shift of the reflected opti- 
cal field. A phase-sensitive measurement of the reflected 
light provides therefore a measurement of the force. 

These optomechanical force detectors have a sensitivity 
which is limited by the thermal noise acting on the mirror 
mechanical degrees freedom, as well as by the more fun- 
damental, unavoidable, quantum noise associated with 
the quantum nature of light, i.e., the phase fluctuations 
of the incident laser beam (shot noise) and the radia- 
tion pressure noise, inducing unwanted fluctuations of 
the mirror position. A compromise between these noises 
leads to the so-called standard quantum limit (SQL) for 
the sensitivity of the measurement 0, 3 . Analogous fun- 
damental limitations affect also other similar detection 
devices, such as nano- and micro-electromechanical sys- 
tems, which are also extensively studied for the realiza- 
tion of ultra-sensitive detection devices p| such as force 
detection on the atto-Newton level and mass detection 
on the zepto-gram level Q . 

Many proposals for the detection of weak forces in- 
volves high-finesse optical cavities with a movable mir- 
ror, in which the phase sensitivity is proportional to the 
cavity finesse J8J. However, recently Ref. j^] has pro- 
posed a new optomechanical detection scheme involving 
a single highly reflecting mirror, shined by an intense 
highly monochromatic laser pulse. A vibrational mode 
of the mirror induces two sidebands of the incident field, 
the Stokes and anti-Stokes sideband. This effect was 
recently observed in a micro-mechanical resonator as a 
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consequence of the radiation pressure force acting on it 
|lfj| . Under appropriate conditions on the duration of the 
laser pulse, the two sideband modes show significant two- 
mode squeezing, i.e., they are strongly entangled [Tlj . 
In particular the difference between the two amplitude 
quadratures and the sum of the phase quadratures of the 
sideband modes can be highly squeezed 113 . and this 
reduced noise properties are used in Ref. [!j to achieve 
high-sensitive detection of a force acting on the mirror. 
However Ref. 9] considered only partially the effect of 
the thermal environment of the mechanical mode. In 
fact, Ref. considered the limiting case of a laser pulse 
duration much shorter than the mechanical relaxation 
time and neglected all the dynamical effects of damping 
and thermal noise. Here we drop this assumption and we 
take into account the effects of the thermal environment 
acting on the mechanical mode, by adopting a quantum 
Langevin equation treatment [13. We shall see that, as 
expected, damping and thermal noise have a detrimen- 
tal effect on the force detection sensitivity, but that one 
can still go below the SQL at achievable values of me- 
chanical damping and temperatures, provided that the 
two sideband modes are appropriately entangled at the 
input. 

The outline of the paper is as follows. In Sec. II we il- 
lustrate the model describing the force detection scheme, 
while in Sec. Ill we define and evaluate the minimum 
detectable force. In Sec. IV we consider experimentally 
achievable parameters and compare the performance of 
the scheme with the SQL for the detection of a force Q , 
while Sec. V is for concluding remarks. 



II. THE MODEL 

We consider the system schematically depicted in 
Fig. 2] It consists of a perfectly reflecting mirror shined 
by a pulsed quasi-monochromatic laser at main frequency 
ljq, linearly polarized in the mirror surface and focused 
in such a way as to excite the Gaussian acoustic modes 
of the mirror, in which only a small portion of the mirror 
around its center vibrates. These modes describe elastic 
deformations of the mirror along the direction orthogonal 
to the surface, and are characterized by a small waist w, 
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where 4>o is the angle of incidence of the driving beam, 
p is the power of the incident beam and Av moc ie is its 
bandwidth, while Avdet is the detection bandwidth. 



A. Including damping and thermal noise on the 
vibrational motion 



FIG. 1: Schematics of the opto-mechanical device to detect a 
force F. A vibrational mode of the mirror (probe) related to 
the ladder operator b and oscillating at angular frequency Q 
is excited by the radiation pressure of the incident laser field 
(angular frequency wo)- Light is scattered into the two first 
sideband modes, the anti-Stokes mode at uo + fl (operator ai) 
and the Stokes mode at u>o ~ f2 (operator a%). The force F to 
be detected acts on the mirror along the horizontal direction. 



a large mechanical quality factor Q and a small effective 
mass M. The motion of the mirror is actually deter- 
mined by the excitation of several modes with different 
resonant frequencies. However, a single frequency mode 
can be considered when a bandpass hlter in the detection 
scheme is used 0] and mode-mode coupling is negligi- 
ble. Therefore we will consider a single mechanical mode 
of the mirror only, which can be modeled as an harmonic 
oscillator with mass equal to the effective mass M and 
angular frequency f2, 
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where X b and P b are position and momentum operators 
of the chosen Gaussian vibrational mode, satisfying the 
commutation rule [X b , Pb] = ih. 

As demonstrated in Ref. , in the case of an intense 
classical incident laser field and neglecting fast terms os- 
cillating at ±f2, the interaction between the chosen vibra- 
tional mode and the continuum of electromagnetic modes 
can be written, in the interaction picture (IP) with re- 
spect to the free Hamiltonian of the system, as a simple 
bilinear Hamiltonian involving the vibrational mode and 
the two first optical sideband modes 
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where b is the annihilation operator of the vibrational 
mode given by b = (iPb + MflX b ) / y/2HM£l, ak are the 
annihilation operators of the optical sidebands Stokes 
(k = 1, angular frequency Uq — and anti-Stokes (k = 2, 
angular frequency loq + ^) modes, and the tilded opera- 
tors are those in the IP. The coupling constants in Eq. @ 



The performance of the three-mode optomechanical 
system described by Eq. @ as a detector of a weak 
classical force acting on the mirror has been studied 
in Ref. 0- In this latter paper, dynamical effects of 
the thermal environment of the mirror were ignored, be- 
cause the dynamics have been studied only in the limit 
of interaction times much shorter than the mechanical 
damping time. Here we consider the more realistic situa- 
tion of non-negligible mechanical damping, which, due to 
the fluctuation-dissipation theorem, implies also consid- 
ering the effect of thermal noise on the mirror vibrational 
mode. These effects are described in terms of a quan- 
tum Langevin equation (QLE) for the vibrational mode 
which, in the first Markov approximation, can be written 
as 13] 
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where the brackets {, } represent the anti-commutator, Y 
is a generic Heisenberg-picture operator, 77 is the damping 
coefficient of the mirror and is the stochastic noise 
force, with correlation function Il3lll5l 
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where kg is the Boltzmann constant and T is the equi- 
librium temperature. Applying Eq. © for the mirror 
position and momentum operators, we get 

P b = -MQ 2 X b + i[H inU P b ] - V X b + £(t), 
implying the following evolution equation for b: 
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If we now move to the IP we get 



b = T \Hint , 61 — — — 
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By virtue of Eq. (J2J we have: 



(8) 



This ratio is usually very low when realistic values are 
taken into account (x ^ 6 ~ 10 5 s _1 , - 10 8 s" 1 ) 
[TT|. Then the last term of Eq. © is much smaller than 
the second term, and can be neglected. Moreover, the 
term b + is counter-rotating and since we have already 
neglected fast oscillating terms at the frequency f2, we 
have to neglect it in Eq. JSJ for consistency. In this way 
we arrive at the final quantum Langevin equation for the 
vibrational mode in the IP 
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where we have defined the damping rate 7 and the scaled 
noise bi n (t) as 
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In the limit of large f2 we are considering, the correlation 
functions of this latter noise term becomes simple. In 
fact, using Eq. 10 and the fact that the factor e is 
rapidly oscillating within the timescales of interest, it is 
possible to derive the following correlation properties of 



(bin(t)bi„{t')) 

<M*)&L(*0> 

(bl(t)b m (t')) 



0, 

(l + ri)S(t-t'), 
n5(t-t'), 



(12a) 
(12b) 
(12c) 



where n = i/( e hil / k sT _ ^ j g |- ne mean thermal vibra- 
tional number at the equilibrium temperature T. To 
state it in an equivalent way, in the limit of large Q, 
the properties of the Brownian noise acting on the vibra- 
tional mode become similar to those of the input noise of 
optical systems. 



B. Exact solution of the dynamics 

The Stokes and anti-Stokes modes are not directly 
sensitive to the damping and noise acting on the mir- 
ror. Moreover they do not undergo additional loss mech- 
anisms since they are traveling waves. Therefore, the 
Heisenberg-Langevin equations describing the dynamics 



of the whole system, in the presence of an additional 
constant force acting on the mirror with dimcnsionlcss 
strength /, are 

a\ = xb (13a) 
b = X a\-0a 2 -2-fb + 2v^b m + infe tnt (13b) 

a 2 = el. 



From these we get the equation for b alone 
b(t) + 2 1 b(t) + e 2 b(t) =g(t), 

where 
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After a straightforward calculation the solution for b 
reads: 
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We have assumed ui real, i.e., 7 < 9, which is typically 
satisfied in the experimentally relevant limit of a high-Q 
vibrational mode. Notice that we reobtain the results of 
Ref. in the limit 7 = 0. The exact expressions for 
the optical sidebands annihilation operators are instead 
given by 
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III. FORCE DETECTION SENSITIVITY 

We now consider the real-time detection of the con- 
stant force / applied to the mirror and determine the 
sensitivity of the considered optomechanical system, by 
evaluating the corresponding signal-to-noise ratio. In op- 
tomechanical devices based on radiation pressure effects, 
one typically performs phase-sensitive measurements on 
the reflected beam (the meter) because the force to be 
detected shifts the mechanical probe determining in this 
way a phase-shift of the field. 

As suggested in we consider an appropriate het- 
erodyne detection [17| of the two sideband modes, corre- 
sponding to the measurement of the operator 



Z, 



ip — e ltp cii - e 
— cos ip 



uzi — aJJ + i sin ip (ai + a^j 



where ip is an experimentally adjustable phase. We 
choose ip — 7r and consider in particular the imaginary 
part of Z n , 
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Using Eqs. one gets 

^ = A 1 (t) Y 1 (0) + A 2 (t) Y 2 (0) + B (t) n(0) 
X-9 

+G(t)+ [ dt' [ dsD(t' -s)Y bin (s), (21) 
Jo Jo 

where Y k {0) = [a*(0) - o£(0)]/2i, k = {1,2,6}, Y bin (s) = 
[hn{s)-b\ n {s)]/2i, 
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The signal is given by the absolute value of the mean 
value of the observed quantity, i.e., S — \{Z^(t))\, while 
the noise corresponds to the square root of the variance 
of the same observable, Af = (Zl{t) 2 ) - (Zl(t)) 2 . Since 
we are considering an open system, averaging means tak- 
ing expectation values with respect to the initial state of 
the system and the environment. In the QLE treatment 
this means averaging over the initial state of our optome- 
chanical system and over the noise bi n (t). 

The natural initial state of the optomechanical system 
is the product state p* o t(0) = 1 0) i (0 1 <S> 1 0} 2 (0 1 <S> , where 
the two sideband modes are in the vacuum state and the 
vibrational mode is in thermal equilibrium with mean 
vibrational number n. 
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However, as suggested in the use of nonclassical 

states, and in particular entangled states of the optical 
modes, could improve force detection sensitivity. For this 
reason we consider the following class of pure initial states 
for the two sideband modes, 
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with s G R, that is, a two-mode squeezed state, repro- 
ducing the usual vacuum state initial condition for s = 
and showing entanglement between the two optical side- 
bands whenever s =/= 0. Notice that when s ^ 0, a nonzero 
incident light power is present not only at the carrier fre- 
quency Wo (po), but also at the two sideband frequencies 
wo ± fi (pi, 2), because power is proportional to sinh 2 s 
[l9|. Therefore, if s is sufficiently large, one could have 
non-negligible scattered light at the additional sideband 
frequencies luq ± 2S7 and interference effects at ujq. This 
however happens only at unrealistically large values of 
two-mode squeezing s. Therefore, we consider not too 
large values of s, so that pi, 2 <C po ~ p and neglect 
these additional effects. 

Using the initial conditions 12.'il) and l|24|) , and the fact 
that (bi n (t)} = 0, one gets that the signal can be written 
as: 
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where G(t) is given by Eq. (|22t 0. while the noise is given 
by the square root of the following variance: 

Af=(9- x f {A\(t)(Y 1 ({)) 2 ) + B 2 (t)(Y b (Q) 2 ) 
+A 2 2 (t)(Y 2 (0) 2 ) + 2A 1 (t)A 2 (t)(Y 1 (0)Y 2 (0)) 

dt' [ dsD(t' - s)Y bin (s) ) )}. (26) 




If we compare these results with the corresponding ones 
of Ref. 13, which considered the same force detection 
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scheme, but in the limit 7 — > (which means neglecting 
the dynamical effects of the environment), we see that 
the signal and noise have the same structure, with two 
important differences. First of all, the time-dependent 
coefficients Ai(t), A 2 (i), B(t) 1 G(t) have a modified ex- 
pression due to the nonzero damping rate 7; moreover in 
the present case, the noise has an additional term, corre- 
sponding to the last line of Eq. i|26|l . Using the definition 
of Yb in (s) and the correlation functions of Eqs. (|12|l . one 
gets (Y bin (s)Y bin (s')) = (2n + l)S(s - s')/4, so that the 
explicit expression of this additional noise term is: 
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which is a positive, non-decreasing function of time for 
any positive t. By using Eqs. (|25|I - H27|I and the following 
initial mean values, stemming from Eqs. (|23(1 - (|24|1 . 
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one can obtain the explicit expression of the signal-to- 
noise ratio 1Z = S/V77. In order to have significant 
results 1Z must be greater than a certain confidence level 
a. This parameter is fixed by the experimenter in ac- 
cordance to his trust in the measuring device; for sim- 
plicity we set here a = 1. The sensitivity or minimum 
detectable input of the device is the minimum magnitude 
of the input signal required to produce an output with a 
specified signal-to-noise ratio. It is easy to see that in or- 
der to obtain 1Z = 1 the sensitivity f m in of the apparatus 
of Fig. ^must be at least equal to ^/M /S. This provides 
the following explicit expression of the minimum force 
detectable with the apparatus at issue: 
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Let us note that / and / m i n are dimensionless quanti- 
ties. The scaling factor to pass to a minimum detectable 
force with proper dimensions F can be obtained from 
Eq. Ijl^b ) and the usual definition of the operator b, giv- 
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IV. 



MINIMUM DETECTABLE FORCE AND 
STANDARD QUANTUM LIMIT 



In this section we study the performance of the force 
detection scheme presented here, characterized by the 
minimum detectable force, Eq. (|29l) . In this respect a 
significative benchmark is provided by the so called stan- 
dard quantum limit (SQL), defined in 0]. 

The SQL represents the optimal sensitivity pertaining 
to an ideal quantum harmonic oscillator when used to 
measure a force applied to it. One typically considers 
only the limitations due to quantum uncertainties, i.e., 
those associated with the oscillator's quantum ground 
state, and assumes zero temperature and no damping 
(shot noise limit). Given an effective interaction time r 
between the force and an oscillator of mass M and angu- 
lar frequency f2, the SQL for the detection of a constant 
force is given by Q 
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This means that in principle Fsql can even become zero 
when r tends to infinite. In practice however the inter- 
action time t cannot be too large. In a realistic setup 
the mechanical damping rate 7 is always nonzero and 
this fixes a first upper limit, t <C I/7. Moreover, a 
perfectly constant force is just an idealization; usually 
one has some typical time scale Tp over which the force 
appreciably changes. As a consequence, it is not conve- 
nient to take t ^ Tp because in such a case the mo- 
mentum change induced by the force may average to 
zero; this fixes a further upper bound for r. So while 
in Sec. Ill Bl we have assumed a constant force, this prac- 
tically means that / does not appreciably vary over the 
typical timescale of the system dynamics, which is essen- 
tially determined by _1 . Hence a reasonable choice for 
r is t ~ O -1 , and we shall set 
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in the expression for the SQL, Eq. I|32|l . Since 7 <C O, the 
above choice is consistent with the requirement t <C I/7. 
Furthermore, we show below that this choice results op- 
timal with respect to the final heterodyne measurement 
of the sideband modes. 

In order to compare the sensitivity of Eq. (|29|l , scaled 
with the factor of Eq. (|31|l , with the SQL for the detec- 
tion of a force, Eq. (|32f> . we choose the parameter regime 



illustrated in Table [I] which, even though challenging to 
achieve, is within reach of current technology. This pa- 
rameter choice gives Fsql = 12.2 x 10~ 18 N. 



Parameter 


Value 


27rc/o;o 


600 nm 


f) 


2tt x 10 7 Hz 


P 


50 mW 


M 


5 x 10" 12 Kg 




10 6 Hz 




10 2 Hz 


7 


1 Hz 



TABLE I: Choice of parameters for the detection scheme of 
Fig. EJ 

In Fig. we plot Log 10 (F/FsQL) as a function of the 
interaction time t (i.e., the duration of the driving laser 
pulse), for different values of the damping rate, 7 = 
{0.01,0.1,1} Hz, corresponding to increasingly darker 
grey curves. The temperature and the two-mode squeez- 
ing parameter s are taken to be zero. We see that the 
minimum detectable force is a quasi-periodic function 
with period 2t:/uj ~ 27r/0, presenting a series of min- 
ima at times t = (2k + 1) tt/lj, with k a positive integer. 
This is due to the entanglement dynamically produced 
by the interaction Hamiltonian of Eq. J2J). In fact, the 
minima are obtained when the two reflected sideband 
modes are factorized from the vibrational mode and are 
in a two-mode squeezed state in which the variance of 
the difference of the amplitude quadratures of the two 
sideband modes, as well as the variance of the sum of 
their phase quadratures, are maximally squeezed below 
the shot noise limit Since the measured ob- 

servable, Z\ of Eq. (|20|l . is just the sum of the phase 
quadratures of the two fields, the minima corresponds to 
the minimum noise, yielding the maximum sensitivity of 
the detection scheme. 

As the interaction time increases, the minimum de- 
tectable force at the local minima increases as well, be- 
cause the effect of damping becomes more and more im- 
portant at longer times. As a consequence, the first mini- 
mum at t\ ~ 7r/0, (t\ ~ 15 ms with the parameter values 
of TableHJ) corresponds to the best possible sensitivity at- 
tainable with the device at issue. In typical situations, 
the time of arrival of the slowly varying force to be de- 
tected is unknown (consider for example the case of the 
tidal force of a gravitational wave). In such a case, the 
best detection strategy corresponds to a pulsed regime, 
in which a precise temporal switch fixes a pulse dura- 
tion exactly equal to t\ for the impinging laser, with a 
repetition rate of the order of the damping frequency 7. 

The dashed line in Fig. [^represents the SQL, obtained 
when F/Fsql = 1- It is apparent that for sufficiently low 
values of 7 the sensitivity goes beyond the SQL. However, 
when 7 assumes the more realistic value of 1 Hz (dark- 
est grey in the figure) the entanglement created by the 
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FIG. 2: Plot of the envelope of Log 10 (F/Fsc>L) (rapidly oscil- 
lating at the angular frequency il) versus the interaction time, 
at four different values of damping, 7 = {0.01, 0.1, 1} Hz, cor- 
responding to increasingly darker grey curves. The other pa- 
rameters are given by Table 1, while the mirror temperature 
is T = and there is no initial entanglement between the 
sideband modes, s — 0. The dashed line represents the SQL. 



dynamics is no more sufficient to go below the SQL, in 
none of its minima. However, as shown by Ref. [fj, there 
is a further resource that can be exploited in order to 
beat the SQL, i.e., the two- mode squeezing of the initial 
state of Eq. (|24|) . This factor represents a sort of "static" 
entanglement that can add its effect to that of the dy- 
namically generated entanglement between the sidebands 
and is able to increase significantly the sensitivity of the 
apparatus. 

Hereafter we concentrate on the first minimum of 
the minimum detectable force of Fig. [21 that is we fix 
ti ~ 7r/0. In Fig. [3 the force sensitivity is plotted versus 
the two- mode squeezing factor s, for four different val- 
ues of the damping rate, 7 = {0.01,0.1,1,10} Hz, and 
again at zero temperature. The dashed line is the SQL. 
As expected, the force detection sensitivity worsens for 
increasing damping. As shown by Fig. [2J when s = 
and 7 = 1 Hz the sensitivity is above the SQL. How- 
ever it goes below the SQL in correspondence of a two- 
mode squeezing coefficient s ~ 1.5 (a squeezing of about 
13 dB). Interestingly enough, at fixed 7, the minimum 
detectable force is not a monotonically decreasing func- 
tion of s, but it has a minimum, meaning that for each 
7 there is an optimal two-mode squeezing value maxi- 
mizing the force detection sensitivity. This feature was 
lacking in Ref. [9(, in which the best possible squeezing 
was the highest one, and is a consequence of the inclusion 
of damping and noise in the model. From a physical point 
of view, this means that once that the interaction time 
t\ is fixed, the input entanglement and the dynamically 
generated entanglement interfere in a nontrivial way, so 
that the optimal sensitivity is achieved at a finite value 
of two-mode squeezing s. 

Finally we study the temperature dependence of the 
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FIG. 3: Plot of Log 10 (F/FsQL) versus the two-mode squeez- 
ing factor s, at four different values of damping, 7 = 
{0.01, 0.1, 1, 10} Hz, corresponding to increasingly darker grey 
curves. The dashed line represents the SQL. The other pa- 
rameters are given by Table 1, while the mirror temperature 
is T = 0. 
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FIG. 4: Plot of Log 10 (F/FsQL) versus the two-mode squeez- 
ing factor s, at four different values of temperature, T = 
{0, 0.03, 3, 300} K, corresponding to increasingly darker grey 
curves. The dashed line represents the SQL. The other pa- 
rameters are given by Table 1. 



sensitivity of the detection scheme. In Fig. 01 we show 
the behavior of the minimum detectable force versus 
s, at four different values of the temperature, T = 
{0, 0.03, 3, 300} K and at fixed damping, 7=1 Hz, while 
the other parameters are again those reported in TableQ] 



We see that up to mirror temperatures of the order of 
few Kelvin degrees, a nonzero value of the input two- 
mode squeezing is able to compensate the detrimental 
effects of damping and thermal noise acting on the mirror 
and one can achieve a detection sensitivity better than 
the SQL. This becomes impossible at room temperature, 
where the minimum detectable force becomes larger than 
the SQL, for any value of s. 



V. CONCLUSION 

We have studied in detail the optomechanical scheme 
for the detection of weak forces proposed in Ref . , based 
on the heterodyne measurement of a combination of two 
sideband modes of an intense driving laser, scattered by a 
vibrational mode of a highly reflecting mirror. In partic- 
ular we have considered the dynamical effects of damp- 
ing and thermal noise acting on the mirror vibrational 
mode, which were neglected in Ref. . The dynamics of 
the system is characterized by a bilinear coupling of the 
two optical sidebands with the vibrational mode, which 
is able to generate significant entanglement between the 
two sidebands for an appropriate duration of the driv- 
ing laser pulse |^,[n|. This condition corresponds to the 
highest signal-to-noise ratio for the detection of a slowly 
varying mechanical force acting on the mirror which, for 
extremely low values of the mirror damping and temper- 
ature, can be better than the SQL for the detection of a 
force At more realistic values of damping and tem- 
peratures, the minimum detectable force becomes larger 
than the SQL and the dynamically generated entangle- 
ment is no more able to counteract the effects of damping 
and thermal noise. However, we show that the presence 
of additional entanglement in the input state of the two 
sidebands may improve the performance of the detec- 
tion scheme and one can go significantly below the SQL, 
even in the presence of non-negligible mirror damping 
and not too low temperatures. For example the SQL can 
be beaten by adopting a vibrational mode with a me- 
chanical quality factor fl/^f ~ 10 and at temperatures 
around 3 K. 
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